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Abstract 

We discuss a particular stringy modular cosmology with two axion fields in seven space-time 
dimensions, decomposable as a time and two flat three-spaces. The effective equations of 
motion for the problem are those of the SU(3) Toda molecule, and hence are integrable. We 
write down the solutions, and show that all of them are singular. They can be thought of as a 
generalization of the Pre-Big-Bang cosmology with excited internal degrees of freedom, and 
still suffering from the graceful exit problem. Some of the solutions however show a rather 
unexpected property: some of their spatial sections shrink to a point in spite of winding 
modes wrapped around them. We also comment how more general, anisotropic, solutions, 
with fewer Killing symmetries can be obtained with the help of STU dualities. 
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There has been much interest recently in using string theory to address the problems 
of cosmology, including the issue of the initial cosmological singularity, symbolizing the Big 



Bang in the Standard Cosmological Model [|TJ-|jT0|- Perhaps the best developed approach 
to date is the Pre-Big-Bang cosmology 0, which originally strove to relate the singularity 
problem with inflation. Briefly, the idea rested on the rolling dilaton field, which would nec- 
essarily evolve toward a region of both strong string coupling and high curvatures, resulting 
in a phase of accelerated expansion of the Universe in the string frame, where some (as yet 
unknown) dynamical mechanism would kick in to prevent further increase of the curvature 
and transform the evolution into a more moderate power(ed)-law expansion. This "graceful 
exit" has been shown to be very difficult to attain; we know that the semiclassical corrections 
alone are not capable of producing it 0, O, O. It turns out that in the sub-Planckian 
region of curvatures the influence of the terms leading to singular behavior in the classical 
theory is too strong. Figuratively speaking, the impending singularity does all the "driving" 
at the scales close to the Planck mass. As a consequence, it has become evident that a 
Pre-Big-Bang phase must evolve to the Planckian energies, where we would need a better 
understanding of the effects of quantum gravity to see if a graceful exit can occur there. 
In practical terms, however, this means that the Pre-Big-Bang cosmology still contains a 
Planck phase, much like the Standard Model it was meant to supplant [^4 -[18 . 



On the other hand, it has been suggested that the singularity can be removed, or at 
least redefined, by the winding mode sources in the Early Universe 0. The winding modes 
are expected to contribute significantly to the effective cosmological dynamics at very high 
energies, or equivalently, at small distances. This is because their energy is generically 
inversely proportional to the radius of the compact direction around which they are wound 
up, and hence it grows rapidly as this radius diminishes. In fact, the presence of the winding 
modes should prevent a complete collapse of the compact directions to the zero radius - the 
growing energy of the winding modes should overtake all the other effects and cause the 
scale factor to "bounce" away from zero. Indeed, this is known to happen in some special 
cases, as investigated in |], ||, [[(J. As a matter of fact, the example of reference [1C| is the 
winding mode in disguise: the pseudoscalar axion field of this solution can be dualized to 
its associated three-form, which is then seen to be purely magnetic: it is equal to a constant 
times the volume form of the three-space. We can then formally compactify the three-space 
to a three-torus, and think of the magnetic axion as a membrane wound around this torus- 
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i.e. the winding mode source. The interpretation of the bounce which occurs in the scale 
factor for this solution then becomes clearer: it arises in response to the increased resistance 
of the membrane to being squeezed down to a point. 

Upon observing this fact, one might wonder what might happen if two axions (equiv- 
alently, two axionic charges) are present. Namely, one might be able to find cosmo logical 
analogs to solitonic string solutions, which are generically nonsingular because of a conve- 
nient redefinition of the singularity. Further, the presence of two axionic charges may lead to 
the stabilization of some of the moduli fields, whose singular behavior often leads to Planck- 
scale curvature singularities. Hence it is interesting to study the details of this mechanism, 
in order to shed some light on the small-scale behavior of winding modes. Surprisingly, the 
model we will discuss here turns out to be completely integrable despite its rather compli- 
cated nonlinear structure. We will recognize the equations of motion for the case of two 
axionic charges as those of the £77(3) Toda molecule, and will therefore be able to write 
down a completely general family of classical solutions. The solutions in turn do not possess 
the properties we have hoped for initially: they all have naked singularities, either in the past 
(like the (— ) branch solutions of the chargeless case 0) or in the future (like the (+) branch 
solutions ||) or both. The phase-space dynamics of the model is considerably more com- 
plex than in the cases without any, or with a single axion charge. Depending on the initial 
conditions, there are solutions whose complete evolution takes place in the whole range of 
coupling strengths, similar to the situation without charges. There are also solutions which 
evolve entirely in the region of strong coupling, like in the single charge case (such behavior 
is expected because of the S-duality relationship of string coupling to the chargeless case: 
g' = c 2 g + d 2 g~ 1 , where c, d are real numbers; since g varies on the whole half- line [0, oo), g' is 
bounded below by 2|cd|). Finally, some of the solutions evolve entirely in the region of weak 
coupling, and have a rather unexpected behavior in that their spatial subspaces oscillate 
around zero volume, in spite of a winding mode wrapped around them. The solutions from 
this subfamily will therefore have both past and future singularities. This is in sharp contrast 
to the chargeless case, where the presence of a mode wound around any spatial section of 
nontrivial homotopy ensures that the section must always have finite volume. The effect is 
a consequence of a nontrivial interaction between the modes, given by the Toda potential, 
as a result of which the dynamics is not just a simple extension of the previous results. 
One should keep in mind, however, that our analysis is classical, and that, in principle, the 



2 



quantum corrections could resurrect the bounce expected from the winding modes. We will 
finally comment on the use of string dualities for generalizing the solutions presented in this 
paper, and in particular point out a connection between our solutions and an anisotropic 
Bianchi II family of cosmologies. The novelty of this result is that not only was isotropy lost 
but also the algebra of the basis forms of the solution has been changed. Such changes may 
turn out to be important for the considerations of the cosmological singularity ||14|| . 

We begin by deriving the simplest nontrivial Toda cosmology, the SU(3) case, as a 
solution of the effective action describing the graviton multiplet in seven dimensions. To the 
lowest order in a', 

S = J £x^ge-*{R + (V$) 2 - lflJ, A + 2A 7 } (1) 

The fields $ and H^ u \ = d^B u \ + cyclic permutations are the seven-dimensional dilaton and 
three-form axion fields respectively. The last term, A7, represents the stringy cosmological 
constant, which can arise from the central charge deficit of the target space. We keep it here 
for completeness, although we will ignore it in the most of what follows. This action can be 
thought of as a consistent truncation of any superstring theory, containing only the NS-NS 
fields in addition to the graviton. We will look for the solutions of the form 

ds 7 2 = -n(tfdt 2 + a{tfdx 2 + b{tfdy 2 
H 7 = Pdx 3 + Qdf (2) 

The axionic equations of motion can be written in form notation as d*H = dH = 0, and so 
it is not difficult to verify that the Bianchi identity requires that P and Q be constant. We 
could now write down the equations of motion obtained from the action above, substitute 
the ansatz and look for solutions. We choose to work in the action, however, dimensionally 
reducing the solution by employing the fact that six coordinates are cyclic, and then varying 
the resulting ID action with respect to the dynamical degrees of freedom, as well as the 
gauge parameter n (the lapse function). The action becomes, up to a ID boundary term, 
and after the redefinition of the dilaton according to $7 = <fi + 3 In a + 3 In b, 

s _ ( dte-*{2k n - -(— + - ^ + — + — } (3) 
J l 2 a 6 b 6 n na? nb 2 > 
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Varying this action with respect to a, b, <fi and n, and then choosing the gauge n = exp( 
we find 

p2 e -20 



h 
k 



2 a 6 
Q 2 e- 2 ^ 



2 6 6 

= 3/i 2 + 3A; 2 -0 2 (4) 

3h 2 + 3k 2 = 2 + 2A 7 e~ 2 * - l -e- 2 \^ + ^) 

2 v a 6 6 6 

where /i and are logarithmic time derivatives of a and b, respectively. If the cosmological 
term were absent^, it is not difficult to see that if one, or both, of the charges vanish, 
the equations of motion are readily integrable. If both charges vanish, we can rewrite the 
dynamics in terms of normal coordinates as a set of free fields ||. When only one of 
the charges vanishes the equations are equivalent to a Hamiltonian system consisting of 



a Liouville theory, or the 577(2) Toda model, and several free fields [10]. Moreover, the 



case with a single charge can be obtained from a chargeless solution by an electric-magnetic 
SL(2, R) duality transformation of the reduced 4D action (which in this case is an exact 
symmetry of the classical 4D action because all the £7(1) 1-form fields vanish identically) 
IIJ. On the contrary, this transformation in its conventional form cannot be employed to 



generate solutions with two charges, because when we reduce the theory to 4D, assuming 
that the internal three-torus is carrying the nonvanishing charge to obtain the other charge 
from the space-time axion, we can see that there arise moduli terms which break the SL(2, R) 
symmetry. Let us show this explicitly. As the solution above has six translational Killing 
symmetries (x and y translations), we could reduce the solution down to 4D by eliminating 
three of the cyclic coordinates. In usual cases, after tedious but straightforward algebraic 
manipulations with the reduced action, we would have rewritten it as: 

S = J d i x^gle~^{R A + (V$ 4 ) 2 + ^Tr(ZVm) 2 - ^H^} (5) 

where the subscript represents the 4D equivalents of the original 7D fields, obtained via the 

Kaluza-Klein dimensional reduction. This correctly describes our case with both charges 

vanishing, or the case when only one of the charges is zero, and we consider the subspace it 

1 Solutions for the case A7 ^ can be found by the application of the technique exhibited in [[[9) to the 
solutions to be presented here. 
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would have wrapped around as internal. The cr-model fields m appear after rearranging the 
scalar moduli fields. The correspondence is given by 

»-(£.:£.) '-(! I) 

where g , b and 1 are 3x3 matrices defined by the internal space components of the 
metric and the axion: g = (#ab) and b = (6ab) ■ The axion field strength can be rewritten 

as 

Hauu\ = d/A^ivX + cyclic permutations (7) 

We have assumed here that all the cross terms (1-forms) between the internal space and 4D 
space-time vanish. To see how the 5-duality comes into play, we need to first conformally 
rescale the metric to the 4D Einstein frame, which is given by g^ v = exp(— $4)^4^, and 
then replace the three-form axion field strength in four dimensions by its dual pseudoscalar 
field. The correspondence between them is 



H^ x = e^^g^ uXp V p B (8) 

The hat here denotes the Einstein conformal frame. If we introduce the complex axidilaton 
field A = B + iexp(— $4) we can rewrite the action ([|) as 

8 = J d 4 x^{R, + jf^f + \Tr{lVmf) (9) 

where f denotes complex conjugation. A careful examination of the equations of motion 
derived from this action shows that they are invariant under the following set of transforma- 
tions: 

A -> A' = °-^4 cb-da = l (10) 

a A + b v ' 

where a, b, c, d are all real numbers. The action also remains invariant under ( |10"D because 
the U(l) gauge fields are absent. The transformations (pi]) combine to form the SL(2,R) 
group. This symmetry group is referred to as the electric-magnetic ^-duality. It is clear, 
however, that the application of this symmetry doesn't offer any special advantage when we 
reduce on the subspace which doesn't carry the charge. It is not difficult to see that in this 
case, the solution is S'-selfdual, modulo trivial rescalings (in direct analogy to, for example, 
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dyonic black holes in string gravity). Thus, in order to see if we can use the SL(2, R) group 
to generate a solution with two charges, we have to carry out dimensional reduction along 
a different path. Namely, to be able to mix the axion and the dilaton in four dimensions 
in a nontrivial manner, we should treat the chargeless spatial section as a part of the AD 
space-time, and the charged spatial section as internal. Let y coordinates parametrize the 
internal space. Because the corresponding 2-form axion field formally depends on an internal 
coordinate, say y 3 , via B = Qy 3 dy l dy 2 (note that it is still invariant under translations 
y 3 — > y 3 + since its variation is proportional to a gauge transformation: SB = ( > dy 1 dy 2 = 
d(( > y 1 dy 2 )) we have to carry out the dimensional reduction more carefully. We cannot group 
the internal components of the B field in the 4D cx-model form - the internal charge Q gives 
rise to an effective dilaton potential in 4D instead. After transforming to the Einstein frame 
and replacing the 4D 3-form field with its Hodge dual as in (|S|), we can write the resulting 
action as 

S = J d x^/g 4 {R 4 + + -Tr(/Vm) - -y-^-} (11) 

The last term obviously breaks the SL(2, R) invariance (the cr-model field m actually contains 
only the scale factor b; the reason that we have chosen to keep the action in the above form 
is to underline that only the last term breaks the SL(2, _R)-duality) 0. 

It is then a little surprising to find that this case is still integrable. We can see this 
as follows: first, we manipulate the equations of motion to rewrite the dilaton equation as 
cf) = —h — k, and integrate this twice to obtain cf) = <po + Ct — Inab (</>o and C are the 
two ensuing integration constants). Next, we note that if we define new variables ip and \ 
according to 

fl3 = i GXP( ^ X " 3Ct ~ 0o) b ' = ^ GXP( ^ " 3Ct ~ 0o) (12) 

and we employ the scaled time variable r = Ct, we can rewrite the independent equations 

2 To get a more intuitive picture how this term arises in four dimensions, we could have split the original 
7D 3-form into two pieces, one living in the effective AD internal space and the other in the spacelike sections 
of the AD space-time. The internal space 3-form could then have been Hodge-dualized with ID Levi-Civita 
tensor, to a 4-form, now living entirely in the AD space-time. Upon the subsequent dimensional reduction 
to four dimensions, this term could have been Hodge-dualized again, now with AD Levi-Civita tensor, to a 
cosmological-like term - i.e. precisely the potential we get above. 

3 Note that when C = these transformations break down. Nevertheless, the system is still equivalent to 
an SU (3) Toda molecule, as can be seen by the substitution a 3 = Pexp(3x/2 — 4>o), b 3 — Q exp(3tp/2 — c/)q). 
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of motion as 

x = e*~ 2x $ = e x ~ 2 ^ 

X 2 + ip 2 - xi> = 6 - e^" 2 * - e x ~ 2 ^ (13) 

These equations are now recognized as the Toda SU(3) molecule; they can be identified 
with the motion of a particle in two dimensions, determined by the Hamiltonian H = 
X 2 + ip 2 — icP + exp( , — 2x) + exp(x — 2-0) with the energy E = 6 (because of the chosen time 
normalization). In order to solve this system, we could invoke the theoretical machinery of 



Lie symmetries and the Lax pair construction [20]; suffice it to say that this leads to the 
solution of the form 



e 



*-$> ie - A * T e* = y £B i e XiT (14) 



where the constants A^, Bi and \ are determined by the initial conditions. For our pur- 
poses, it is sufficient to use this as another ansatz, simply substitute it in the equations of 
motion (|T3|) and use the initial conditions exp(x ) = ot\, exp(?/> ) = «2, exp(x )xo = fix and 
exp^o^o = fii to determine the above parameters. 

We will merely sketch here the main points of this computation. Choosing r = as the 
instant to define the initial conditions, we see that combining them with the ansatz ( [Lip 
leads to the identities = «i, Si-Sj = «2, Yl\\Ai = —fix and X^Aj-Bj = fi 2 . It is 

convenient to introduce new parameters ak and bk according to 

a k = (Ai - A fc )(Aj - \k)A k b k = (A, - A fe )(Aj - X k )B k i^j^k (15) 

For now, we will assume that Aj ^ Xj,i ^ j. We will obtain the degenerate case later, 
as a limit of the nondegenerate one. Treating, for example, Aj, a 3 and 6 3 as independent 
parameters we can express the remaining four parameters in terms of them. This split 
however is completely arbitrary and hence permuting the indices we obtain the general 
formulae 

a k - dj = (Xj - X k )(Xiax + fix) b k - bj = (Xj - X k )(Xia 2 - fi 2 ) (16) 

Substituting these identities into the Toda molecule equations of motion gives further con- 
straints, which relate a k 's and b k s: 

3 3 3 

aj bj = i rja^n^^-i X> = o ( 17 ) 
i i i 
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We can now combine the constraints ( |16D and fllTf ) to determine and solving these 
equations yields 

a, = ^-# hj = h°±±h. (18) 

3 AjQJl + A J A-,0! 2 -/3 2 

The final step is to determine the "eigenvalues" Xj. It is straightforward to see that they 
are given as the solutions of the system of equations 

f> = 
i 

3 

n(V*i + #0 = -(a 2 /?i + ai/3 2 ) (19) 
i 

3 

J|(AiQ;2 - (3 2 ) = (® 2 (3i + ai(3 2 ) 



A more transparent notation can be achieved if we introduce the three parameters K n = 
h Z)i(Aj) n , n = 1, . . . , 3 (note that K\ = 0). With the help of these, we can see that the A, 
are the three roots of the secular polynomial 

A 3 - k 2 \ - k 3 = (20) 

What remains to determine is the relationship between the parameters «j and the initial 
conditions, in order to specify their values. This we find from the last two constraints of 
(|19|) . The result is 

K 2 — „ + 9 + 9+9 

af ol\(x 2 a<2 af 

(p 2 ai - /5ia 2 ) + -o 9 ( 21 ) 



«3 



9 9 Vr^-^ -L r'i-^ZJ i 9 9 



We recognize the first of these two equations as k 2 = H = 6 - i.e. this is nothing else but 
the Hamiltonian constraint in a different guise. The second equation can be rewritten as 
K 3 = jokO* ~ X) + X ex P(x — 2 , 0) — ^ exp(-^ — 2%); this is also an integral of motion as can 
be verified by differentiating it with respect to time, and using the equations of motion. It 
is a consequence of a hidden symmetry of the Toda system, and effectively the reason why 
the SU (3) Toda molecule is integrable. 
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We can therefore summarize the solution as follows: 



e 



A 1 A fc a 2 - & 
* " £ A ^ x ^ + eAfcT ( 22 ) 



^ A(A fc ) A fc a 2 - /3 2 



A(A fc ) = E II (A* -A,) 

Note that a simultaneous change of sign of the roots A& and a permution of the fields x an d 
"0 amounts to time reversal. This will be useful below. At this point we could substitute this 
solution back into the ansatz for the graviton multiplet (|2p. A more transparent representa- 
tion, in our opinion, can be obtained by a firsthand qualitative analysis of (0), which would 
highlight the physically important aspects of an Universe described by (|2|) and (|22|). One 
should note that whereas the solution fl22|) is valid for all real r, the definition (fj) of x an d 
■0 implies that the times for which the scale factors vanish may be singularities, and we have 
to examine them carefully. Let us first eliminate unphysical parameters from the problem. 
It is easy to see that the effect of the initial conditions a\ and a 2 on the bulk dynamics 
is ignorable; by the coordinate transformations dx k — > dx k / P^^Ohj 2 , dy k — ► dy k /Q 1 / 3 ^ 
we can set the charges P,Q equal to «^ 3 ^ 2 , a^ 2 , respectively. Therefore, the effect of the 
initial conditions for x and ip is only to set a possible deficit angle if the spatial sections are 
three-torii. In what follows we will therefore set ai — a 2 = 1- 

Now we can continue with the discussion of the solutions. There are three subfamilies 
described by (|22|), depending on the roots A& of the secular polynomial. The possibilities 
are i) roots all real and different; ii) roots all real, two of them degenerate; Hi) one root 
real, two complex (conjugates of each other). The roots cannot all be imaginary because 
of the positivity of the Hamiltonian constraint. One could expect now that the functions 
exp(x) and exp(^) are always positive, in analogy to the case with a single charge. This is 
indeed true in cases i) and ii); the Universes described by these two subclasses of solutions 
are still singular, undergoing first a long period of contraction until the volumes become 
small enough for the winding mode source to overturn the contraction of at least one (but 
not necessarily both) of the scale factors into a stage of unbounded expansion. Depending 
on the sign of the integration constant C, these solutions have either a past (C < 0, (— ) 
branch) or a future singularity (C > 0, (+) branch). The case Hi), however, has a completely 
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different singularity structure, where the spatial sections always shrink to a point in spite of 
the charges they carry This is without an analogue in the simplest models of Pre-Big-Bang 
cosmology. 

We now present the analysis case-by-case. In case i) we will first demonstrate that the 
functions exp(x) and exp(^) are positive, as claimed above. When all the roots are real 
and different, we can order them as a decreasing sequence; so suppose Ai > A2 > A3. For 
simplicity, we can assume k 3 > 0. This implies that Ai > and A 2 ,A 3 < (because 
K 3 = n^fc)- The other case, k 3 < 0, can be reduced to this one by a time inversion (i.e. 
the sign change of A^'s and the interchange of x an d 4>)- So we have Ai > > A2 > A3. If 
P1 + P2 7^ 0, the Eqs. ( pp tell us that a/- 7^ aj for k 7^ j. For suppose the opposite; if any two 
a fc 's were equal for different fc's, we would have obtained the identity (/?i+/3 2 )(Afe — Xj) = for 
the pair of indices k,j under scrutiny; this is a contradiction as the roots are nondegenerate 
and we have assumed (3\ + p 2 7^ 0. Further, aia 2 a 3 = —1 tells us that two of the a^'s are 
positive and one is negative. We now need to show that ai,a 3 > and a 2 < 0; this will 
guarantee that A^s and B^s, and hence exp(x) and exp(ip), are always positive. We will 
do this by deriving inequalities for the linear functions A& + Pi and A& — p 2 . There are four 
special cases to look at: a) Pi,p 2 > b) /3i,/3 2 > c) j3\ > 0, p 2 < with f3\ + (3 2 > 
and d) /?! < 0,/3 2 > with (3% + /3 2 > 0. The options with the opposite sign of (3% + (3 2 
are equivalent to the last two cases after the interchange of \ and ip (which implies the 
interchange of /3i and [3 2 and hence the sign flip of (3\ + (3 2 ). So we find: 

a) by our assumptions, Ai + /3i, (3 2 — X 2 , and (3 2 — A3 are all positive; from n&=i(Afc+/3i) = 
]lfe=i(/^2 — Afc) = —fix — (3 2 < and A 2 > A 3 we see that A 2 + f3\ must also be positive, while 
A3 + Pi and p 2 — Ai must be negative; this follows by the negativity of the product of these 
functions, and the fact that A3 + (3\ > and A 2 > A3 would imply A 2 + f3\ > 0, contradicting 
the assumption that their product is negative; using these inequalities in the definitions of 
cik we immediately see that a 1; 03 > and a 2 > 0; in conjunction with A(Ai), A(As) > and 
A(A 2 ) < 0, we find Ak, Bk > as claimed. 

b) the products satisfy I\l =1 (Xk + Pi) = Yll=i(^ ~ ^fc) = ~ A — P2 > 0; by assumptions, 
A 2 + Pi, A 3 + Pi and p 2 — Ai, must be negative; this and A 2 > A 3 imply that Ai + Pi and 
P 2 — A 3 are positive, whereas p 2 — X 2 is negative; for suppose that P 2 — A 2 is positive - this 
would imply that p 2 — A3 > 0, contradicting the assumption that the product of these two 
expressions is negative; again, this means that a 1; 03 > and a 2 > and so Ak, B^ > 0. 
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c) this time, ft > 0, ft < and ft + ft > 0; thus IlLi(** + ft) = IlLi(#2 - A fc ) = 
—ft — ft < 0; since ft and ft + ft are the same as in a) above, we have Ai + ft > 0, 
A 2 + ft > and A3 + ft < 0; by the assumption about signs, ft — Ai < 0; therefore we see 
that either ft — A 2 > 0, ft — A 3 > or ft — A 2 < 0, ft — A 3 < 0; the latter situation would 



imply ai, a 2 > and a 3 < 0; but by a 3 — a 2 = (A 2 — A 3 )(Ai + ft) (see Eq. ([16])) we would 
get a 3 > a 2 , contradicting the latter case; hence it must be ft — A 2 > 0, ft — A 3 > and so 
ai, a 3 > and a 2 > again yielding A k , B k > 0. 

d) finally, ft < 0, ft > and ft + ft > 0; so njU = E[fc=i(^a-A fc ) = -ft-ft < 
0; since now ft and ft + ft are the same as in b), we get ft — Ai < 0, ft — A 2 < and 
ft — A 3 > 0; the signs imply that Ai + ft > 0; the remaining options are A 2 + ft < 0, 
A3 + ft < or A 2 + ft > 0, A 3 + ft > 0; the latter case is again contradictory, because it 
would lead to a\, a 2 > and a 3 < whereas by a 3 — a 2 = (A 2 — A 3 )(Ai + ft) we would get 
a 3 > a 2 ; hence it must be A 2 + ft < 0, A 3 + ft < and so a\, a 3 > and a 2 > 0, once more 
leading to A k , B k > 0. 

To completely cover the case of real nondegenerate roots, we need to review the situation 
when ft + ft = 0. Fortunately this is simple enough to be evaluated directly The equality 
ft = —ft = (3 and the Hamiltonian constraint set ft = 4/3. Moreover, the condition 
Ilfc = i(Afc + ft = tells us that one of the roots must be equal to —ft So we factor the 
cubic secular polynomial, to find the roots: — ft ±\^E + (3/2 with (3 = ±2/v3. Ordered 
roots corresponding to (3 = — 2 /a/3 are \/5 + l/y/3, — 2/y / 3 and — \/5 + l/\/3 (the /? of 
the opposite sign is equivalent to this one by a permutation of our degrees of freedom, as 
in previous cases). Clearly a\ = a 3 = 1, but because A 2 + j3 = 0, the formula for a 2 is 
ill-defined, leading to an expression of the form oc 0/0. A more careful calculation should 
employ a 2 — a 3 = (A3 — A 2 )(Ai + ft = — 2 to give a 2 = —1. Lastly, using the values for A& 
and the formulae for A k , B k , we find A\ — B\ — l/2y / 5(v / 5 + V / 3), A 2 = B 2 = 1/2 and 
A 3 = B 3 = 1/2^5(^5 - a/3) - i.e. A k , B k > 0. 

In sum, we see that when the roots are real and nondegenerate, the functions exp(x) 
and exp(^) remain positive, each being a linear combination of three positive terms. This is 
reminiscent of the case of a single charge and a constant internal space radius [ lOj] , whereby 
we might be tempted to conjecture that the spatial sections of this solution can never shrink 
to a point. Intuitively this would agree with the assertion that a winding mode causes an 
initially contracting scale factor of a subspace around which the mode wraps to bounce into 
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expansion. The analogy is misleading, however, because the scale factors are not identical 
with the solution of the Toda equations and we need to carry out a more careful inspection 
of the solutions in order to see what happens. 

First of all, the solutions are still singular, as can be seen from the inspection of the Ricci 
scalar. A straightforward computation yields 

R = 6e 2 * (h+'k + (j)(h + k) + 2h 2 + 2k 2 + 3hk) (23) 

Let us look at this expression for the case of three real nondegenerate roots when r — > ±00. 
As before, Ai > > A2 > A3. So in the limits we have 

e x -> -^-e- A3T -> * N e Air as r^oo 

A(A 3 ) aiA(Ai) 

e x -> -r^-e- AlT -> — x e A3T as r -> -00 (24) 

A(Ai) a 3 A(A 3 ) 

The expression for in terms of the scale factors a and b allows us to rewrite the Ricci scalar 
as follows: 

R = 6C*e^e^(*±i + X 2 + ^~^ _ x _ ^ (25 ) 

The limits ( p4|) imply x — > —A3, ^ — >• Ai when r — ► 00, and \ — > — Ai, ^ ^ A3 when 
r — > —00, with x, ip — > in both cases. The Ricci curvature therefore behaves as R oc 
exp((6 =F (Ai — A 3 ))r) as r — > ±00, while remaining finite for all other times. Using the 
constraints for Ki and k 2 (Eqs. (^9|)-(pT|)), we can show that 6 =F (Ai — A3) > 6 — \[2A > 0. 
Thus, as r — ► 00, the curvature diverges, and as r — ► —00, the curvature vanishes^. Further, 
we see from the expressions for the Hubble parameters that in the limit r — > — 00 both scale 
factors are decreasing - i.e. the spatial sections are collapsing. This follows from the limits 
of the logarithmic derivatives of scale factors, which can be evaluated from (|1^) and (|24|): 
h — > — (Ai + 2)/2 and k — > — ( | A3 1 + 2)/2. The limit r — >• 00 is more ambiguous; to study it, 
let us substitute the solutions and the choice of the lapse function back into the ansatz for 
the metric. In the limit r — > 00 we find (recall that A3 < 0) 

where K T , K x , K y are numerical constants irrelevant for the discussion that follows. As 
we have commented above, 6 — Ai + A 3 > 0, and so the appropriate comoving time is 
4 The same is true for the special case C = mentioned in footnote 2. 



12 



T = exp((Ai — A 3 — 6)r/2), and in the limit r — > 00, T — > 0. In terms of it, the metric is 

rfs 2 ^ _^ rfT 2 + ir^rff 2 + K yl ^-dy 2 (27) 

where ni = 2(|A 3 | - 2)/(6 - Ai - |A 3 |), n 2 = 2(Ai - 2)/(6 - Ai - |A 3 |). Clearly, the signs of 
these two numbers determine the qualitative behavior of the solutions: for rik positive, the 
scale factor expands superexponentially, and for n k negative it shrinks to zero according to a 
simple power law. This kind of behavior is of course familiar from the simplest Pre-Big-Bang 
scenario; what is novel is the combination of these possibilities for the given values of charges. 
The allowed combinations are simultaneous pole expansion of both spatial sections (as, for 
example when Ai = —A3 = v^6) and the contraction of one spatial section accompanied with 
the pole expansion of the other. The simultaneous contraction of both spatial sections to a 
point is impossible; this can be seen as follows. The Hamiltonian constraint, n 2 — 6 tells us 
that the roots A fc live on a sphere of radius y/l2 : J2^l = 12. Using Ai > > A 2 > A 3 and 
assuming simultaneous contraction of both scale factors (and hence Ai < 2, |A 3 | < 2), we get 
J]A^<3x4 = 12- which is impossible. The last ingredient of the phase-space description of 
evolution is the behavior of the coupling constant. In the limits r — > ±00, the coupling of the 
theory exp ($ 7 ) = exp(0)a 3 6 3 behaves asexp($ 7 ) = exp(x+^> — 3t) — > exp((±(Ai — A 3 )— 3)r). 
From this we see that in the limit r — > — 00 the coupling diverges for all values of A&, while 
for r — > 00 the coupling diverges when Ai + |A 3 | > 3 (all the cases of simultaneous expansion 
of both scale factors, and some of the cases with one scale factor increasing and the other 
decreasing), and vanishes when Ai + |A 3 | < 3 (the remaining cases when one spatial section 
asymptotically superinflates and the other shrinks to a point). 

With these details, a coherent picture emerges at last: the phase-space of the problem 
can be divided into three different sectors, distinguished by the behavior of the solution 
near the singularity. In one sector, as the singularity is approached, the evolution is in the 
strong coupling regime with both scale factors growing beyond limit. In the other sector, 
the coupling is still strong, but instead of simultaneous expansion, one of the spatial sections 
shrinks to a point. In the last sector of the phase space again one of the scale factors grows 
and the other decreases, but the system evolves towards the weak coupling regime. All 
the solutions are smoothly connected to their other limit, where they approach an almost 
flat space where both spatial sections contract simultaneously, and the string coupling is 
in the strong coupling regime. By time inversion (reversing the sign of C), we see that 
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both directions of evolution are admissible: we either start from one of the singular regions 
and evolve towards an asymptotically flat space ((— ) branch solutions) or begin from a 
strongly coupled asymptotically flat space and progress towards one of the singular sectors 
((+) branch solutions). This generalizes the simplest versions of the Pre-Big-Bang scenario. 
The difference is that, depending on the initial conditions, the solutions studied so far are 
qualitatively similar to both the chargeless (solutions evolving throughout the allowed range 
of variation of the string coupling) and the single charge case (solutions whose dynamics is 
confined to the strong coupling regime). The bottom line is that the singularity cannot be 
removed at the classical level from any solution characterized by a set of real nondegenerate 
roots. 

The case with real degenerate roots is qualitatively similar to the previous case. We will 
here prove only that the functions exp(x) and exp(^) never vanish for any finite value of r. 
This is again the essential property of the solution - the rest follows straightforwardly. We 
will view this limiting situation of the previous one. In order to quantify this, let 



us introduce a parameter e = A2 — A3, rewrite the solution (22) in terms of it, and take the 
limit e — > + . To do this, note that the roots A& can be rewritten as Ai = A — e, A2 = e — A/2 
and A3 = —A/2 (A is the nondegenerate root of the secular equation). Then we can write 

3 



e 



fc Ji a fc (e)A fc (e) 

where A^. are at most quadratic polynomials in e, as can be seen from their definition in 
([Z2|). Also, by their definition, the coefficients at are linear fractions in e. Hence the solution, 
being a finite linear combination of a finite product of functions analytic in e must also be 
analytic in it; so we can take the limit e — > + , which is unique and reproduces the case 
of two degenerate roots. It is then easy to verify that the result of taking the limit can be 
written as 

Xr ( , Xr/2 

A x (0) 



cti(0)Ai(0) 

where p(r) and q(r) are some linear functions of r, p(r) monotonically decreasing and q{r) 
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monotonically increasing. To see that exp(x) and exp('0) never vanish, it is sufficient to 
prove that the lines p{r) and q(r) never intersect the curves exp(— 3A/2r) and exp(— 3A/2r). 
The flow of the curves is consistent with this assertion: the lines generally stay below the 
exponential curves. We still need to exclude the possibility for any intersections, or in fact 
even osculations of the lines with the exponentials. The following qualitative argument does 
precisely that. Consider the case when e is very small but not zero. Suppose then that 
any of the two functions of the solution in the limit e = 0, say exp(x), vanishes for some 
r. Analyticity of exp(%(e)) ensures that we can choose a sufficiently small negative e such 
that exp(x(e)) < 0. But this would contradict the fact that exp(x) and exp(ip) are always 
positive for real nondegenerate roots, proven above. Hence we conclude that the same must 
be true when two of the roots are degenerate. The rest of the analysis for this case can now 
proceed along the same lines followed in the case of nondegenerate roots, leading to similar 
conclusions. 

The last case we need to consider is when two roots are complex (conjugates). This case 
is dramatically different from the two previous situations, as we will now show. We can write 
down the roots as Ai = A, A2 = —A/2 + ifi and A3 = —A/2 — ifi. The solutions then become 



Introducing the trigonometric notation and defining pi = |a 3 /A 3 |, p 2 = l/|a 3 A 3 |, 8\ = 
Arg a 3 /A 3 and 82 = — Arg a 3 A 3 , we can rewrite these equations as 



It should be obvious from the these expressions that for a sufficiently large (small) r, exp(x) 
(exp(?/>)) must vanish. The reason lies in the oscillatory nature of the cosines, which for some 
t can, and will, turn negative enough to cancel the exponentials. However, all the zeros of 
the expressions in fl3~T|) are simple. This can be verified by looking at the function f(z) = 
rjexp(z) + ( cos(z) and computing the period around one of the zeros: n = {l/2m) § df / f = 
1. An even simpler way to see this is to plot the cosine and the exponential against each 




(30) 




(31) 
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other and compare the tangents at the intersection points - the fact that they are different 
guarantees that the roots are simple. As a corollary, the first derivative of \ must have a 
simple pole with unit residue at the location of the root of exp(x). Also, we see from ( |3lD 
that the other function from the set of solutions can neither vanish nor diverge here. We 
can now reexamine the Ricci curvature for this case, and with the help of the equations of 
motion rewrite it as 

R = 6C 2 e^ 0+6r + — + V*-^ + io\) - e-x-^x ~ e"*"^) (32) 

In the limit when r approaches the location of the root of exp(x), the term of the leading 
order of divergence is exp(— 3%); it blows up as 1/ (t — r ) 3 , as our simple analysis of the zeros 
and poles combined with power counting shows. Hence there is a curvature singularity at 
the location of any of the roots of exp(x), and similarly, of exp (?/>). In response, we conclude 
that the case of two complex (conjugate) roots is pathological: all solutions have both a 
past and a future singularity, where one of the scale factors must vanish. We also note that 
this subfamily of solutions is divided into two sectors: in one, exp(ip) is always nonzero, 
while exp(x) vanishes twice in the lifetime of the Universe, and vice versa. These two groups 
of solutions are divided by a single solution straddling r = 0, where the past singularity 
corresponds to the first root of exp(x), and the future singularity to the first root of exp(^) 
(or vice- versa, depending on the sign of of C). Finally, the string coupling is again given 
by exp(<3> 7 ) = exp(x + tjj — 3r). When the solution approaches the singularity, the coupling 
vanishes, and at all other times it remains finite. Thus, the solutions from this subfamily 
evolve between two consecutive zeros of scale factors, always remaining in the weak coupling 
region. 

If we take any of the solutions above as a starting point and reduce it to 4D such that 
the resulting metric is of the Robertson- Walker form, we will see that all those modular 
cosmologies still suffer from singularities. The only difference may be in deforming the phase 
space in certain ways - for example, redefinitions of effective 4D couplings may change the 
strength of couplings in the phase space of the lower-dimensional theory. Still, no curvature 
singularity can be removed in such a way, as can be noted from the fact that the reduction 
does not affect the 4D space-time metric, and that in all the solutions discussed above 
the curvature singularity is encoded democratically in both scale factors. Nevertheless, the 
reduction to fewer dimensions may permit us to utilize the duality symmetries for extending 
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the solutions in nontrivial ways. We will merely quote here the possibility of reducing the 
solution on one of the spatial directions to QD, and using either the U duality between 
the IIA and the IIB theory, or the SL(2,R) NS-NS/RR-RR duality in IIB theory as a 
solution-generating technique. This would of course give rise to Toda models carried by 



forms of different rank (e.g. a 2-form and a 1-form, as mentioned very recently in ||23|| ). 
We will not delve into a detailed investigation of such solutions here. Instead we will point 
an interesting relationship between our solutions, reduced to AD, and a class of modular 
Bianchi II cosmologies. If we look at the metric and the axion of the AD version of our 
ansatz, ignoring the internal degrees of freedom, 

ds^ 2 = —n(t) 2 dt 2 + a(t) 2 dx 2 
B 4 = Px l dx 2 A dx 3 (33) 
$4 = $7 -3 In b 

we see that the translations along x 2 and x 3 are manifest symmetries of the solution (the 
translation in the x l direction is also a symmetry since the change of the axion 2-form is just 
a gauge transformation, as discussed on page 6). So we can transform this solution with the 
simplest of all T-dualities, the R—>l/R map plj . We can choose to dualize with respect to 
x 2 , in which case the map is 

, _ 1 , _ B 2ll , _ 92^ 

i/22 — i/2/x — JD 2a — 

#22 922 <?22 

9^u = 9n» (34) 

^22 

r,! c 92^B 2u — B 2fJi g 2u 

B - = B "" S 

The dual solution can be written as 

ds' A 2 = -n(t) 2 dt 2 + a(t) 2 ((dx 1 ) 2 + (dx 3 ) 2 ) + -^(dx 2 + Px l dx 3 ) 2 

a[t) 

B' 4 = $^ = $ 4 -21na (35) 

Without engaging in a detailed analysis of this solution, let us merely point out that it belongs 
to a different class of homogeneous, but not isotropic cosmologies: the Bianchi II Universe, 
defined by the following algebra of the spatial basis 1-forms e 1 = dx 1 , e 2 = dx 2 + Px 1 dx 3 
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and e 3 = dx 3 : 

de 1 = de 3 = de 2 = Pe 1 A e 3 (36) 

This change of the algebra of basis 1-forms reflects the change of the topology of space-time 
under the duality fl34|) . We find this relationship rather peculiar. While both solutions, in 
their classical form, are singular, it is known that the a! expansion resummation of certain 
anisotropic cosmological models can lead to the removal of the singularity present at the 



classical level, and thus to the Planck-scale resolution of the graceful exit problem [14, 17, 18 



A similar argument for homogeneous and isotropic models is still lacking, in part because 
of the absence of a conformal field-theoretic description needed for such an investigation. 
Perhaps such an approach may be aided by a relationship of the type considered here. 
The T-duality, and perhaps the other dualities of string theory, may provide a natural way 



to introduce such anisotropies |22| into the consideration of apparently homogeneous and 
isotropic models, and so allow for a resolution of the singularity even in this case. 

In conclusion, we have shown how a nontrivial SU (3) Toda molecule dynamics arises in 
string cosmology from a gravity-mediated interaction of two winding modes. The model is 
exactly integrable, and leads to a broad family of solutions which are all singular, either in 
the future or in the past, or both. Thus all solutions suffer from the graceful exit problem, if 
viewed as a realization of the Pre-Big-Bang cosmology. In addition, there are solutions with 
string winding modes wrapped around spatial subspaces, where the subspaces can shrink 
to zero volume. This aberration from the behavior anticipated earlier comes about as a 
consequence of the gravity-mediated interaction between the modes. We have also pointed 
out a duality relationship of our solutions with a Bianchi II anisotropic model, and have 
suggested that while the model is classically singular, the anisotropies may lead to its Planck- 
scale regularization, after the a' expansion is resummed. The relationship we have pointed 
out could also be viewed as an indication that anisotropies in string cosmology may not have 
as drastic consequences for the evolution of the Universe as they do in General Relativity 
- for if they can be removed altogether by a duality map, believed to relate equivalent 
physical pictures, they need not be a physically relevant notion. The clarification of the 
relationship between the singularity, anisotropies and duality therefore seems to warrant 
further investigation. 



18 



Note added in proof: While this report was being prepared, there appeared an article 
by A. Lukas, B. Ovrut and D. Waldram, ||23|| , which overlaps with a portion of the work 
presented here. 
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